The a-axis optical conductivity of detwinned ortho-II YBa 2 Cu 3 O6.50 
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The a-axis optical properties of a detwinned single crystal of YBa2Cu3O6.50 in the ortho II phase 
(Ortho II Y123, T c = 59 K) were determined from reflectance data over a wide frequency range (70 
- 42 000 cm -1 ) for nine temperature values between 28 and 295 K. Above 200 K the spectra are 
dominated by a broad background of scattering that extends to 1 eV. Below 200 K a shoulder in 
the reflectance appears and signals the onset of scattering at 400 cm" 1 . In this temperature range 
we also observe a peak in the optical conductivity at 177 cm -1 . Below 59 K, the superconducting 
transition temperature, the spectra change dramatically with the appearance of the superconducting 
condensate. Its spectral weight is consistent, to within experimental error, with the Ferrell-Glover- 
Tinkham (FGT) sum rule. We also compare our data with magnetic neutron scattering on samples 
from the same source that show a strong resonance at 31 meV. We find that the scattering rates 
can be modeled as the combined effect of the neutron resonance and a bosonic background in the 
presence of a density of states with a pseudogap. The model shows that the decreasing amplitude 
of the neutron resonance with temperature is compensated for by an increasing of the bosonic 
background yielding a net temperature independent scattering rate at high frequencies. This is in 
agreement with the experiments. 

PACS numbers: 74.25.Gz, 74.62. Dh, 74.72. Hs 



The complete phase diagram of the high temperature 
superconducting (HTSC) cuprates is still under intense 
debate. The normal state, particularly in the underdoped 
region, is dominated by a variety of not-well-understood 
cross-over phenomena that may either be precursors to 
superconductivity or competing states. These include 
the pseudogap Q, the magnetic resonance SSB, the 
anomalous Nernst effect [fj , stripe order H, Jala E3 
and possible superconducting fluctuations |lll Il2j ] . The 
situation is further complicated by the presence of dis- 
order and the practical considerations that lead to a sit- 
uation where a given cuprate is not investigated with 
all the available experimental techniques. Ideally, one 
would like to have a system where disorder is minimized 
and several experimental techniques can be used with the 
same crystals. As a step in that direction we present de- 
tailed a-axis optical data on the highly ordered ortho-II 
phase of Yf^CuaOg.so (Ortho II Y123) and compare 
these data with recent results from magnetic neutron 
scattering and microwave spectroscopy on crystals from 
the same source. 

An important motivation for a comparison between 
transport properties and the magnetic neutron reso- 
nance comes from the observation that the carrier life 
time, as measured by infrared spectroscopy, is domi- 
nated by a bosonic mode HI Q H [H E3 whose 
frequency and intensity, as a function of temperature 
and doping level, tracks the inelastic magnetic reso- 
nance at 41 meV with in-plane momentum transfer of 
(ir, 7r) 0, 0| . The magnetic resonance has also been in- 
voked to explain other self-energy effects such as the kink 
in the dispersion of angle-resolved photo emission spec- 
tra (ARPES) HlEllSilim and as a hump-peak struc- 



ture in tunneling spectra |22j . These effects have also 
been attributed to the electron-phonon interaction [23j| . 
While the bosonic excitation may not be the fundamen- 
tal engine of superconductivity [la, yj| , nevertheless it is 
important to map out the regions in the phase diagram 
where it can be found and correlate the various experi- 
ments that yield evidence of its presence. 

The YBa2Cu307_a; (Y123) material has been one of 
the most thoroughly studied of all the HTSC systems 
but, like most cuprates, it suffers from disorder associ- 
ated with the charge reservoir layer, in the case of Y123, 
disordered oxygen chains. However, at an oxygen dop- 
ing level of x — 0.50, an ordered ortho-II phase occurs 
with alternating full and empty CuO chains, doubling 
the unit cell along the a-direction, and yielding a well 
ordered stoichiometric compound [251 |2fj . A further ad- 
vantage of this system is the availability of very large 
crystals, suitable for neutron scattering. A disadvantage 
of the Y123 system is that the crystals are not easily 
vacuum-cleaved and as a result, surface sensitive probes 
such as angle resolved photoemission or scanning tunnel- 
ing microscopy have been used less with Y123. In con- 
trast Bi2Sr2CaCu208 (Bi-2212) cleaves easily and unlike 
Y123, can also be overdoped. Optical spectroscopy has 
the advantage of working on both systems equally well 
and offers a bridge between the ARPES surface sensi- 
tive experiments and the large-volume probe of neutron 
scattering. 

The paper is organized as follows. We start with a 
brief introduction to the experimental method followed 
by a presentation of the raw reflectance data and an anal- 
ysis where we use the extended Drude model to extract 
the various optical constants, in particular the real and 
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imaginary parts of the scattering rate. To improve the 
accuracy of the data, we first remove known features in 
the conductivity spectra such as the transverse optical 
phonons and the delta-function response of the super- 
conducting condensate. The measured scattering rates 
are then compared with a theoretical model. We next 
discuss the overall results by comparing them with data 
from other experiments on Ortho II Y123, in particular 
recent neutron dataji3| on samples from the same source. 

EXPERIMENTAL METHOD AND RESULTS. 

The detwinned ortho-II YBa2Cu30g.5o sample used in 
this study was grown by a flux method using BaZrC>3 
crucibles HsJ, annealed under pure oxygen gas flow at 
760 °C, and detwinned at 300 °C by a pply ing uni-axial 
stress of 100 atm along the a-direction [251] . The dimen- 
sion of the sample is 1.5 x 1 x 0.2 x mm 3 . The nearly 
normal-incident reflectance of the sample was measured 
over a wide frequency range (70 - 42000 cm -1 ) at nine 
temperatures between 28 K and 295 K with a Bruker IFS 
66v/S Fourier transform spectrometer with linearly po- 
larized light. A polished stainless steel mirror was used 
as an intermediate reference to correct for instrumental 
drifts with time and temperature. An in situ evaporated 
gold film on the sample was the absolute reflectance ref- 
erence [2^| . The reflectance of the gold films was in turn 
calibrated with a polished stainless steel sample where 
we relied on the Drude theory and the dc resistivity as 
the ultimate reference. An advantage of this technique 
is that it corrects for geometrical effects of an irregular 
surface. The in situ gold evaporation technique gives 
an absolute accuracy of the reflectance of better than 
±0.5 %. For the reflectance data above 14 000 cm -1 we 
used aluminum instead of gold as the coating material. 
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TABLE I: The fitting parameters of Ortho II Y123 reflectance 
at three temperatures, T = 28 K, 100 K, and 295 K. The unit 
of LUpj, LUj, and -yj is cm -1 . An extra band for T = 28 K is 
the super fluid band at zero frequency. The extra band for 
T — 100 K shows additional spectral weight at low frequency 
at low temperature (see Fig. 0. The peak SI splits at 100 K 
into two components. 
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FIG. 1: The a-axis reflectance of Ortho II Y123 at nine tem- 
peratures. There is a broad band of absorption giving rise to 
a region of negative curvature of the reflectance between 350 
and 600 cm -1 . At high temperature the curvature is posi- 
tive at all frequencies. The inset shows the high frequency 
reflectance. 
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FIG. 2: The low frequency reflectance is shown as thin lines 
and an oscillator fit is shown as heavy lines. The data are fit- 
ted with a Drude band and a series of oscillators, one of which 
is a strong electronic mode marked SI. In the superconduct- 
ing state, an oscillator at zero frequency is added to represent 
the superconducting condensate. The oscillator parameters 
are shown in Table I. The open circles near zero frequency are 
reflectance data at 28 K from a microwave measurement l33l. 



Figure^shows the raw reflectance data for nine differ- 
ent temperatures in the frequency region from 75 cm -1 
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to 3000 cm" 1 . The inset extends the data up to 42 000 
cm" 1 . A strong temperature dependence is observed in 
the low frequency range but becomes weaker at higher 
frequencies. Two specific features show temperature de- 
pendent amplitudes: a broad shoulder at » 400 cm" 1 
and a less obvious depression in reflectance at 175 cm" 1 . 
In addition, signatures of well-known transverse optical 
phonons can be seen as sharp minima in the low tempera- 
ture spectra At low temperature the overall 
curvature of the reflectance is negative at low frequen- 
cies changing to a positive curvature at high frequency, 
which gives rise to an inflection point at ~ 700 cm" 1 . 
The inset shows the high frequency reflectance with a 
prominent plasma edge at 13 000 cm -1 and a relatively 
weak temperature dependence above 25 000 cm -1 . 

Figure [21 shows the lowest frequency range on an ex- 
panded scale. In addition to the transverse phonons, we 
observe a new feature at 175 cm" 1 , that appears be- 
low 171 K, and is denoted SI. The thin lines are the 
measured reflectance curves and the solid thick lines are 
oscillator fits where the phonons have been modeled as 
Drude-Lorentz oscillators. The electronic background is 
represented as a broad Drude peak in the normal state 
and a two-fluid model with an additional delta function 
at the origin in the superconducting state. The param- 
eters used in the fit are shown in Table I. From the fit 
we get a superconducting condensate density of 72 x 
10 6 cm -2 expressed as the square of a plasma frequency 
in cm" 1 . This corresponds to a penetration depth of 
1900 A. In comparison, the recent Gd Zero-field Electron 
Spin Resonance (ESR) measurements on samples from 
the same source yield a penetration depth of 2000 A and 
are shown as open circles at very low frequency at 28 
K |3^ . Our estimated overall error in the absolute value 
of the reflectance of ± 0.5 % translates to a ± 300 A error 
in penetration depth. We conclude that our data are in 
agreement with the ESR results, consistent with errors 
associated with the measurement of absolute reflectance. 
If we assume an overall monotonic variation of the fitted 
curves with temperature, we can estimate, from the devi- 
ation from this monotonicity, a relative error of 0.2 % for 
the temperature dependence of the reflectance. Finally, 
from the noise level of the spectra, estimated to be 0.05 
% at 300 cm" 1 and rising to 0.3 % at 100 cm -1 , we can 
set an upper limit of 250 cm -1 to the plasma frequency 
of any phonon features we could have missed above 200 
cm" 1 . This limit would be higher at lower frequencies. 

The optical conductivity and the other optical con- 
stants were determined from the measured a-axis re- 
flectance by Kramers- Kronig analysis j3^|, for which ex- 
trapolations to lo — ► and 00 must be supplied. For 
lu — ► 0, the reflectance was extrapolated by assuming a 
Hagen-Rubens frequency dependence in the normal state, 
(1 — R) oc a; 1 / 2 , and below T c an (1 — R) cx uj 4 extrapola- 
tion was used. We used the oscillator fits shown in Fig. [21 
to investigate the sources of error arising from the low 
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FIG. 3: The upper panel shows the optical conductivity where 
several prominent phonon peaks have been identified as PI, 
P2 etc. In the lower panel the phonon peaks (PI to P4) have 
been removed. A complex set of absorption bands can be seen 
at low temperature dominated by a broad peak at 177 cm" 1 
designated as SI. 



frequency extrapolations. The reflectance was extended 
to high-frequency (between 40 000 and 350 000 cm" 1 ) us- 
ing data from Romberg et al. |3^ • Free-electron behavior 
(i? oc oj~ a ) was assumed to hold at higher frequencies. 

The top panel of figure [3 displays the optical conduc- 
tivity for the nine temperatures, eight in the normal state 
and one in the superconducting state. The optical con- 
ductivity can be written as <j(u>) — —iu>(e(u>) — £h)/47t 
where cr is the dielectric constant at a high frequency 
(~ 2 eV). At all temperatures we observe four sharp 
phonon modes out of the six infrared active phonon 
modes [3(3, E01 expected for this polarization of the 
incident radiation. The phonons have been labelled PI 
- P4 in the figure. We fitted the optical conductivity of 
these lines with a Drude-Lorentz model and found the 
line intensities to be temperature independent and of a 
magnitude expected for transverse optic (TO) phonons. 
Table I shows the parameters of the model. 

The bottom panel of figure[2|shows the optical conduc- 
tivity without the four phonon modes obtained through 
the subtraction of the fitted Drude-Lorentz conductiv- 
ities from the measured conductivity. Several features 
stand out. First, there is a prominent onset of conduc- 
tivity that appears at low temperature around 400 cm" 1 . 
This feature is common to all cuprate superconductors 
and was identified early on as an onset of scattering from 
a bosonic mode0, Q, |3(|. Secondly, a strong peak 
at w 180 cm -1 , designated SI, grows as the tempera- 
ture is lowered. At lower frequencies, below 150 cm" 1 , 
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FIG. 4: The partial spectral weight up to 2000 cm" shown 
as dashed and solid curves. The spectral weight increases 
monotonically with temperature and frequency except for the 
28 K curve (heavy solid curve) in the superconducting state 
where there is a loss of spectral weight to the condensate 
and the curve is sharply lowered. The curve marked N s — 
N e f f {67 K) — Nef f (28-RT ) is an estimate of the spectral weight 
of the superconducting condensate. The curve marked N' s is 
a plot of the condensate density obtained independently from 
the imaginary part of c(w). 

there appears to be an additional strong absorption band, 
but any structure here must be interpreted with cau- 
tion since they are derived from reflectance data close to 
unity which are subject to large systematic errors. Fi- 
nally, at higher frequencies, there is a broad continuous 
background absorption that extends up to the plasma fre- 
quency and has been attributed to the influence of strong 
correlations |37U38| . 

We define the effective number of carriers per cop- 
per atom in terms of the partial sum rule: N e ff{uS) = 
2 "^2 U Jq+ <J\{uj')duj' where m is the free electron mass 
and Vcu is the volume per copper atom, 57.7 A 3 . Fig- 
ure 0] shows N e ff{oj) calculated from the optical conduc- 
tivity. In the frequency region shown N e ff(u>) increase 
uniformly with frequency and temperature in the normal 
state. 

In the superconducting state at 28 K, shown in Fig- 
ure 0] as a heavy solid curve, there is a dramatic loss 
of spectral weight due to the formation of the super- 
conducting condensate. The effective number of carri- 
ers per copper atom in the condensate can be estimated 
from the partial sum rule by subtracting the normal state 
curve at 67 K from the superconducting curve at 28 K, 
also shown in the figure. The resulting difference curve 
rises rapidly at low frequency and saturates at a value of 
N s = 0.0264 above 1000 cm" 1 . To find a better value of 



FIG. 5: The partial spectral weight integrated up to various 
frequencies as a function of temperature. Below 16 000 cm" 
there is an increase in spectral weight as the temperature is 
lowered signaling a line narrowing on this frequency scale. Be- 
low T c there is strong loss of spectral weight to the supercon- 
ducting condensate. There is no evidence of any precursors 
to superconductivity at 67 K. In the inset we show N e ff{io) 
at 295 K. 

N s we need to estimate the normal state curve at 28 K. 
We do this by extrapolating the temperature dependence 
in the normal state to 28 K (see Fig.lSJ. With the cor- 
rected value of N e f / at 28 K we get a more accurate value 
of N s = 0.0276. The condensate density can also be esti- 
mated from the imaginary part of the conductivity pro- 
vided the contribution from the residual conductivity is 
not included i.e. N' s = ^^aj[a 2 {uj) - a 2 r(u)] = 0.0283 
where a 2r is the residual quasiparticle conductivity in the 
superconducting state. Figure H shows both N s and N' s 
and it is clear that the two curves (when corrected) ap- 
proach one another closely and the two methods differ 
by 2.5 % well within our 5 % estimated experimental er- 
ror. There have been several reports of ab-plane spectral 
weight changes on enterin g th e superconducting state, 
in other high T c cuprates |H HJ 01 at the 1 to 2 % 
level, but it should be noted that systematic errors in 
the absolute value of reflectance have a large influence 
on the magnitude of both N s and N' s , although the er- 
ror of the difference is smaller. We estimate that our 
error in N s and N' s is 30 % but less than 5 % in the 
difference N' s — N s . Thus, at this point, we are unable 
to conclude from our data that there is evidence for any 
added or missing low frequency spectral weight when the 
superconducting state forms in the Ortho II Y123. 

Figurc0shows the temperature dependence of the par- 
tial spectral weight at various frequencies referred to the 
spectral weight at 300 K. We note that below 16 000 
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cm -1 the spectral weight increases as the temperature is 
lowered while above this frequency there is a decrease. 
At 16 000 cm -1 the partial spectral weight is temper- 
ature independent. Since this frequency is close to our 
estimate of the limit of the free carrier conductivity, we 
conclude that there is a transfer of spectral weight to low 
frequencies as the temperature is lowered within the free 
carrier band on two frequency scales: a lower scale of the 
order of 500 cm -1 shown by the spreading of the curves 
in Figure 0] and a higher one of the order of 5000 cm -1 
that is responsible for the temperature dependence below 
10 000 cm -1 in Figure^ We also note here that the tem- 
perature dependence of the partial spectral weight shows 
no sign of any precursors to superconductivity at 67 K 
which is 8 K above the bulk superconducting transition 
temperature. 

The extended Drude model offers a detailed view of the 
charge carrier scattering spectrum and its contribution to 
the effective mass ^i} . In this picture the scattering rate 
in the Drude expression is allowed to have a frequency 
dependence: 
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FIG. 6: The optical scattering rate obtained from the ex- 
tended Drude model. Two onsets, denoted by arrows domi- 
nate the scattering. 



CT(W ' T) 1 4ttuj+[u J \{uj,T) + i/t(uj,T)} 

= 4 l - (1) 

where lo v is the plasma frequency 1/t(oj,T) is the scat- 
tering rate and A(w)+1 = m*(w)/m, m*(uj) is an effective 
mass and m the bare mass. We also introduce the optical 
self energy S°p = E° p + iS° p , where -2£° p = u\(w,T) 
and — 2E2 P = V T - The optical self energy is, apart from 
a cos{9) — 1 factor, where 9 is a scattering angle, an 
average over the Fermi surface of the quasiparticle self- 
energy [l(| [HI 00] as measured by ARPES. Contrary 
to expectations, in optimally doped Bi-2212 where both 
optical and ARPES data exist, the self-energies derived 
from the two s pec troscopies |45j are surprisingly similar 
to one another |15l Il6l Il8| . 

Figure shows the scattering rate calculated 
from the extended Drude model where 1/t(lu,T) = 

-2E° p (w, T) = ^Re(^j T y) for Ortho II Y123. We also 
evaluate the real part of the optical self-energy, which is 
given by: -2£°>, T) = u,A(w, T) = -glm(^) -w 
and is shown in Figure |S1 For the calculation of the 
optical self-energy and the scattering rate we need a 
value for the plasma frequency which includes spectral 
weight up to the interband transitions. Using a proce- 
dure adopted in a previous study we find a plasma 
frequency w p =15 110 cm -1 . 

In Fig. we see an overall increase in scattering both 
with temperature and frequency. The frequency depen- 
dence is monotonic with two thresholds, shown with ar- 
rows, where the scattering rate undergoes a step-like in- 
crease, a prominent low frequency one at 400 cm -1 and 




FIG. 7: The optical effective scattering rate. Given by the ra- 
tio of the real and imaginary parts of the conductivity times 
the frequency, this quantity has a linear frequency dependence 
that extends to very high frequencies and a temperature de- 
pendence that does not vary much with frequency. 

a weaker, higher frequency one at 850 cm -1 . The vari- 
ation of the scattering rate with temperature is roughly 
linear at low frequency, but becomes much weaker at high 
frequency above the high frequency threshold. The scat- 
tering rate is larger than the frequency at all tempera- 
tures except at 28 K in the superconducting state, i.e. 
Ti/t > huj. 
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FIG. 8: The optical self-energy. The imaginary part of the 
scattering rate is plotted. A broad peak is seen at 295 K (the 
lowest curve). As the temperature is lowered this peak splits 
into two components with a sharp peak appearing below 175 
K as shown by the curvature change at this temperature. The 
inset shows a curvature analysis of the self-energy where the 
slope of the self-energy in the two regions, shown in the main 
panel as two dashed lines, have been plotted, circles, high 
frequency and diamonds, low frequency. 



Fig.[7|shows the effective optical scattering rate denned 



as: 



1 



(2) 



One advantage of this quantity is that to calculate it the 
plasma frequency is not needed. This quantity is also 
more linear with frequency than the optical scattering 
rate |47j . 

In Fig. [S] we show the real part of the optical self- 
energy. It is related by Kramers-Kronig transformations 
to the scattering rate shown in Fig. El The separation of 
the two components of scattering is less obvious in this 
plot but in analogy with the Bi-2212 system where we 
have investigated the detailed doping dependent data on 
this quantity 0, we find that the bosonic mode gives 
rise to a peak on top of a broader background. With 
increased underdoping in Bi-2212 the peak becomes tri- 
angular and harder to resolve from the background and 
is very similar in appearance to the Ortho II Y123 data 
shown here. We can attempt to resolve the peak by fo- 
cusing on the break in slope that occurs between 1000 
cm -1 and 1400 cm -1 in the temperature range from 28 
K to 171 K. Above 200 K the break in slope cannot be 
resolved. The inset to Fig. OD shows a plot of the slopes 
of two straight dashed lines that have been fitted to the 
experimental data above and below the break frequency. 
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FIG. 9: The Eliashberg function, a 2 F(uu) » W(cu), which is 
calculated by using the Eq. [3] We observe two peaks: one is 
temperature dependent at lower frequency and the other is 
temperature independent at higher frequency. In the super- 
conducting state there is a large shift of the low energy peak 
and a dramatic overshoot of the scattering rate giving rise to 
a negative contribution to the Eliashberg function. 



We find that a plot of the two slopes shows that they 
cross at 155 K, i.e. the second derivative of the optical 
self-energy goes through zero at this temperature. This 
is shown in the inset to Fig. |5| 

Recent calculations [2l] of the self-energy of the charge 
carriers interacting with collective spin excitations sug- 
gest that the self-energy acquires an S-shaped frequency 
dependence in the presence of a sharp mode, whereas in 
the absence of such a mode the curve has a monotoni- 
cally negative second derivative. From these observations 
we conclude, as a first approximation, that the bosonic 
mode is confined to temperatures below 155 K. 

For another estimate of the strength of the interac- 
tion of the charge carriers with the sharp mode we used 
the procedure introduced by Marsiglio et al. fl5t l36l lisf 
where the bosonic spectral function is derived from the 
second derivative of the optical scattering rate. This 
function can be written as follows |3fj, |42j, |49| : 



W{u) 



1 d 2 

2ltduJ 2 



(3) 



and W(u>) « a 2 F(uj) at zero temperature in the nor- 
mal state, where a is a coupling constant, and F(u>) is a 
bosonic density of states. The results are shown in Fig.El 
The normal state shows two broad peaks; a prominent 
one at w 350 cm -1 and a much weaker one at w 800 
cm -1 . The higher frequency peak does not change with 
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FIG. 10: Area under the SI mode as a function of temper- 
ature. We also show the area under the neutron resonance 
peak H3] at Q = (1/2, 3/2, 2.2) and 33.1 meV, solid triangles 
from 25 to 43 meV (open triangles) as a function of temper- 
ature. 



FIG. 11: DC resistivity extracted from the optical conduc- 
tivity by extrapolating the optical conductivity to zero fre- 
quency. The temperature dependence shows the familiar lin- 
ear variation with an intercept on the temperature axis of « 
90 K. 



temperature while the lower frequency one grows mono- 
tonically as temperature decreases. Another interpreta- 
tion is in terms of a temperature dependent peak and a 
temperature independent background, four times lower 
in amplitude. In the superconducting state there is a 
dramatic shift of the lower peak to higher frequency and 
the development of a region of negative spectral function 
W(ui) between 700 and 1000 cm -1 . This behavior has 
been predicted by Abanov et al. [isf . 

Fig. El shows the area under the peak SI as a function 
of temperature obtained from a fit of a Lorentzian func- 
tion to the conductivity peak. We see that the intensity 
of this feature decreases linearly with temperature going 
to zero at rj 200 K. On the same graph we have plot- 
ted the area under the neutron resonance and it is clear 
that the two phenomena have parallel temperature de- 
pendencies. We note that the spectral weight of the SI 
peak, 8000 cm -2 represents only 0.006 % of the total free 
carrier spectral weight, 1.25xl0 8 cm' 2 . 

We note that to obtain a finite frequency absorption 
peak such as the peak SI, we need to break transla- 
tional symmetry of the charge carriers. We can rule out 
phonon absorption which would be inconsistent with the 
large spectral weight of the peak SI and the fact the 
spectral weight is temperature dependent as shown in 
Fig El There are a number of possible candidate mech- 
anisms but we would like to focus on the static charge- 
density patt erns seen by scanning probe microscopy in 
Bi-2212 |50j,|51|. Models for optical absorption by strip- 
like patterns have been described by Benfatto and Morais 



Smith 10]. The model describes the optical absorption 
by charges confined by a static pinning potential with 
two parameters, the strength of the potential and the 
modulation amplitude. 

Finally, we show the calculated dc resistivity extracted 
from the extrapolated zero frequency limit of the opti- 
cal conductivity. Figure ^] shows the resulting dc re- 
sistivity for the eight normal state temperatures. The 
resistivity shows the familiar linear temperature depen- 
dence with an intercept on the temperature axis of 90 
K, typical of undcrdopcd samples of high T c cuprates. 
A more detailed comparison shows some discrepancies 
with the best dc transport data J^. First, our absolute 
value of the resistivity at 300 K is 0.7 (mSlcm) -1 while 
Ando et al. find, at the same doping level, a resistivity 
of 1.0 (mficm) -1 . The second difference is the value of 
the resistivity in the low temperature plateau region just 
above T c where Ando et al. measure a residual resistivity 
of about 200 (^Slcrn)" 1 whereas our extrapolation yields 
less than 50 (fiilcm)^ 1 . We think this difference may 
originate in the more ordered chain layer in our samples 
which were annealed whereas the sample of Ando et al. 
were quenched and presumably more disordered. The 
lower absolute value of the 300 K resistivity is harder to 
understand since we don't expect inelastic scattering, a 
mechanism responsible for the high temperature resistiv- 
ity, to depend on disorder. One often suspects contact 
geometry in dc transport but as Ando et al. have shown, 
errors from those sources amount to at most a few per- 
cent in their data. Where the two data sets are in good 
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agreement is the value of ~ 90 K for the intercept on the 
temperature axis. 



DISCUSSION 

We now turn to a general interpretation of our data 
in terms of models proposed for the electrodynamics of 
cuprates focusing our attention on the low frequency 
spectral region below the free carrier plasma edge « 10 
000 cm -1 . It is clear from even the most superficial view 
of the data that the optical properties shown here are 
not those of a set of free carrier interacting weakly with 
impurities through elastic scattering or with lattice vibra- 
tions through inelastic scattering. Two approaches have 
been used to account for these deviations, the one and 
the two component models of conductivity. In the one 
component model any deviation from the simple Drude 
formula are attributed to a single set of charge carriers 
interacting with a variety of bosonic excitations while in 
the two component model a low frequency Drude com- 
ponent is separated by curve fitting and any remaining 
absorption is attributed to a mid-infrared band. 

There is ample evidence that many conducting oxides 
do have a separate well-defined mid-infrared band with a 
peak frequency that moves to lower frequencies with dop- 
ing |5^| . This is also true of the one-layer cuprates, in 
particular L^-zSr^CuC^ |54|. In some materials, when 
the peak reaches zero frequency, striking changes occur 
in the transport properties, for example, Bai-^K^BiOa 
becomes superconducting at this doping level |5Jj|. Re- 
cent evidence from thermal conductivity suggests that 
L^-zSr^CuC^ becomes an insulator below a doping level 
of 0.06 where superconductivity also sets in whereas Y123 
remains in the conducting pseudogap state well below 
the doping level of the superconductivity boundary. 
These observations suggest that the one-component to 
two-component transition occurs at different doping lev- 
els for these two materials and that the Y123 system re- 
mains metallic to the lowest doping levels. We conclude 
that, at least at the level of doping of ortho-II Y123, we 
are justified in using the one-component model. 

The striking result of the one-component model is 
the strong frequency dependent scattering rate 1/t(u>, T) 
shown in Fig|B| There is a nearly perfect linear variation 
of scattering at high frequency with a strong threshold at 
400 cm" 1 at low temperatures that gradually washes with 
increasing temperature. We will interpret our spectra 
within the one-component model in terms of a scattering 
rate that is averaged over the Brillouin zone and is the 
sum of two frequency dependent terms, a bosonic mode 
that dominates at low temperature but whose spectral 
weight weakens as the temperature rises and a feature- 
less temperature independent background that extends 
to high frequencies. 

Several theoretical models of strongly correlated sys- 



tem have been evoked to explain the strong linearly ris- 
ing frequency dependent scattering rate. Earlier mod- 
els included the Mar gina l Fermi Liquid^^i the Nested 
Fermi Liquid models |58j |. the Nearly Antiferro mag entic 
Liquid |33,|5^|, and the Luttinger Liquid model [33 • All 
models predict a continuous linear rise of scattering gov- 
erned by an energy scale of the order of w 0.5 eV. The 
large energy scale rules out a simple phonon mechanism 
for the broad background scattering which would predict 
a flattening of the scattering rate beyond the maximum 
phonon frequency ^3 • The strongly temperature depen- 
dent threshold has generally been interpreted in terms of 
coupling to a bosonic mode 11311141, in particular the the 
41 meV magnetic resonance |l5j, |la, LL2J ■ 



Bosonic Mode Analysis of the Optical Scattering 
Rate 



Here we will adopt the approach of Schachinger et 
al. and treat both the mode and the continuum scat- 
tering on an equal basis assuming that both originate from 
coupling of the charge carriers to bosonic fluctuations. 

To study the temperature dependence of the coupling 
of bosonic modes to charge carriers we used the follow- 
ing expression to model the scattering rate within the 
extended Drude formalism: 



+ OO 



dn a 2 F(n) dz[N(z-n) + N(-z + n)} 



[» B (fi) + 1 - f(z - n)} [f(z - u) - f(z + u)] 

(4) 



where a 2 F(£l) is the bosonic spectral function [flcj . 
N(z) is the normalized density of state of the quasipar- 
ticles, n B (Q) = l/(e m - 1), f(z) = l/(e /32 + 1) are the 
boson and fermion occupation numbers, respectively, and 
/3=l/(k B T). 

Eq.0] represents a finite temperature generalization of 
the T = expression j^H > 



2tt 



dfta 2 F(Q) 



dz-[N(z) 



■N(-z)] 



(5) 

Eq. H|is derived using the method proposed by Shulga 
et al. |62l | and its derivation will be given elsewhere |63l | . 
Both Eq. 0] and are suitable for the case when the 
quasiparticle density of states, N(z), cannot be taken as 
constant in the vicinity of the Fermi levels, e.g. in the 
pseudogap state. 

The density of state N(z) was modelled with a pseu- 
dogap with a quadratic gap function: 



N(z) = [N(0) + (l-N(0))^]6(A-\z\)+6(\z\-A) (6) 
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where 1 — iV(O) is the depth of the gap, 6(z) is the Heav- 
iside step function, and A = 350 cm -1 is the frequency 
width of the gap based on spectroscopic data for the pseu- 
dogap from tunneling |r3^|. and c-axis infrared conduc- 
tivity 65] . An example for a depth of the gap of 0.50 is 
shown in Fig. 

For the bosonic spectral function a 2 F(£l) we used the 
sum of two functions, a peak and a background. 



a 2 F(n) 
PK(Q) 

BG(fi) 



PK(fi) +BG(f7) 
A 



{n-n PK y/[2(d/2.35y] 



27r(d/2.35) 



n 2 



(7) 
(8) 



(9) 



where PK(f2) is a Gaussian peak and BG(f2) is the back- 
ground which we have modelled on the spin fluctuation 
spectrum of Millis, Monien, and Pines (MMP) H^. A 
is the area under the Gaussian peak, d is the full width 
at half maximum (FWHM), and flpK is its center fre- 
quency. We fixed the parameters of the Gaussian peak 
to values shown in Tabled based on the inelastic neutron 
scattering results of Stock et al. [6(|. I s is the intensity 
of the MMP background and Slo is the frequency of the 
background at maximum. The complete bosonic spec- 
tral function a 2 F(£l) used is shown in Fig. It has 
two adjustable parameters, the amplitude of the Gaus- 
sian peak A and depth of the pseudogap 1 — N(0). All 
the other parameters have been determined from other 
experiments. 

We used least squares to fit our scattering rate data 
shown in Fig. © to Eq. g| The amplitude of the MMP 
background and its center frequency were determined by 
fitting the data at 295 K including only the MMP back- 
ground in a 2 F(tt) with a very shallow gap (see Table ITT|) 
in the Fermion density of states N(z). We fixed the back- 
ground parameters for all other lower temperatures to 
their 295 K values. For further fits at lower tempera- 
tures, only two free parameters were used: the depth of 
the gap in the density of states and the area under the 
resonance peak, A. The calculated 1/t(u>) spectra are 
compared with the measured data in Fig ll5l The dimen- 
sionless coupling constant or mass enhancement factor A 
is defined as: 



A 



Jo 



a 2 F(n) 

n ' 



(10) 



The contributions to Xpk and XbGi from the peak and 
the background as well as other parameters of the model 
are shown in the Table IIT1 

One result of our fit is the temperature dependence 
of the depth of the pseudogap shown in the lower panel 
of Fig. El The depth, we find, corresponds to a 41 % 
depression of the density of states at low temperature 
decreasing gradually to 17 % at room temperature. The 




320 640 960 1280 1600 1920 

Frequency (cm' 1 ) 

FIG. 12: The bosonic spectral function a 2 F(Q.) obtained from 
the least square fits to the experimental data. In the inset we 
show determined by neutron scattering from Fig. 14 

in Ref. 66]. 




Temperature (K) 

FIG. 13: The upper panel shows the quasiparticle density of 
state with a gap, N(z). The lower panel shows the depth of 
the pseudogap for the least square fits. The other parameters 
of the fit are shown in Table [n] 



depth of the pseudogap in the c-axis optical conductiv- 
ity is more pronounced 65] but it is known that c-axis 
transport is weighted more heavily in the antinodal di- 
rection where the pseudogap is deeper whereas ab-plane 
transport is more evenly distributed in momentum space. 

Next, we compare our spectral function a 2 F(fl) in 
Fig. E] with the magnetic susceptibility x"(w) deter- 
mined by neutron scattering at 6 K (< T c ) by Stock et 
al. (see inset in Fig. 112ft . The two sets of curves are very 
similar although it should be pointed out that our data 
are in the normal states and neutron data are in the su- 
perconducting state and that the width of our mode has 
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been set equal to that of the neutron resonance mode. 

Comparing the relative amplitudes of the peak and the 
background, we find a substantial difference with the neu- 
tron data. For the ratio of the background to the peak 
amplitude our a 2 F(fl) at 67 K gives ~ 0.17 while the 
neutron \"{ui) gives ~ 0.27 in the normal state (we esti- 
mated the ratio in the normal state from the temperature 
dependent intensity of the neutron mode 0| ) . The rela- 
tive amplitude of the peak in our spectra is 63 % stronger 
than what is seen by neutron scattering. In the absence 
of a more detailed calculation of the spectral function, 
this discrepancy is not surprising since the coupling of 
the charge carriers to spin fluctuations involves integra- 
tions over the Fermi surface and we are reporting only 
on weighted averages. On the other hand our data leaves 
open the possibility that the high frequency channel of 
conductivity, that is active in producing our background, 
is not well described by the MMP spin fluctuations. In 
addition to the relative amplitudes we can also compare 
the areas under the peak and background spectral func- 
tions. There is a sum rule for the area under x"(w) which 
is related to our optical function a 2 F(il) by a coupling 
constant, g 2 . Neutron scattering finds that the area un- 
der the spin resonance is about 3 % of the total area. In 
our case at 67 K this fraction is about 25 %, considerably 
larger. It can be argued that the coupling to the reso- 
nance g 2 could be larger than its background value by a 
factor of ~ y/8. This is reasonable since the resonance is 
around (tt,tt) where the susceptibility is also expected to 
peak. 
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FIG. 14: The mode contribution (upper panel) and the MMP 
back ground contribution (lower panel) on the scattering rate 
based on Eq. 0] with the parameter shown in Table [H] and 
pseudogap information of Fig 1131 for the normal states. We 
note at high frequencies the opposite temperature dependen- 
cies of the mode and the background contributions. 



temp. 


peak 


MMP 


depth 


(K) 


A Qpk d Xpk 


I s Qo Xbg 


1 - N(Q) 


67 


255 248 80 2.20 


354 320 3.30 


0.413 


100 


235 248 80 1.93 


354 320 3.30 


0.413 


147 


170 248 80 1.40 


354 320 3.30 


0.386 


200 


105 248 80 0.86 


354 320 3.30 


0.411 


244 


62 248 80 0.51 


354 320 3.30 


0.366 


295 


248 80 0.00 


354 320 3.30 


0.172 



TABLE II: The parameters of the bosonic mode analysis at 
six representative temperatures, T = 67 K, 100 K, 147 K, 
200 K, 244 K, and 295 K. The peak is a Gaussian function 
(see Eq. 8) and MMP is the background (see Eq. 9). The 
quantities Xpk and Xbg are the coupling constants for the 
peak and the MMP background, respectively. The depth, 1- 
N(0), is the depth of the gap at Fermi level in the density of 
state. All the frequencies are measured in cm -1 . 

In Fig. E| we compare three quantities: the total area 
under the magnetic susceptibility in the 25 to 43 meV 
energy range from Stock et al. |2JJ, the area under the 
peak at 350 cm -1 in the W(ti) fa a 2 F(n) obtained from 
the second derivative of l/r(w), and the area under the 
peak at 248 cm -1 in the a 2 F(£l) from our fit. We note 
that the two areas from the two different procedures for 
deriving a 2 F(fl) show almost identical temperature de- 



pendencies. The area under the peak obtained from the 
magnetic susceptibility shows a temperature trend simi- 
lar to those of the optical data. 



DISCUSSION AND CONCLUSIONS. 
The two channels of scattering of the charge carriers 

In analyzing the optical properties of Ortho II Y123 we 
have adopted the one-component model of conductivity, 
attributing any deviations from the simple Drudc formula 
to inelastic processes. We find that the high frequency 
processes are well described by the MMP spectrum of 
spin fluctuations with a temperature dependence that 
originates in the bosonic occupation numbers, but a tem- 
perature independent spectral function, while at lower 
frequencies a strong temperature dependence sets in and 
the conductivity is well described by the extended Drude 
model with a temperature dependent spectral function. 

The two methods we used to estimate the shape and 
temperature dependence of the bosonic spectrum gave 
very similar results. The main contribution to the spec- 
tral weight came from a prominent peak at 350 cm -1 
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FIG. 15: The optical scattering rates (solid lines) and our 
fits (symbols) based on Eq.[I]with the bosonic spectral func- 
tion shown in Fig. 1121 and pseudogap information of Fig. 1131 
for the normal states. We note that the negative tempera- 
ture dependence of the mode contribution combined with the 
almost equal and opposite contribution of the background re- 
sults in a nearly temperature independent scattering rate at 
high frequencies, in complete agreement with the experimen- 
tal results. 
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which was strongly temperature dependent in strength. 
From this temperature dependence we identified the 
mode with the neutron resonance which has a frequency 
of 248 cm -1 . The 100 cm" 1 discrepancy we attributed to 
the effect of the pseudogap in the density of states which, 
as calculations suggest, has the effect of shifting upward 
the frequency of any bosonic mode. To fit the data, the 
strength of the mode had to be temperature dependent 
and the electronic density of states had to have a gap 
whose depth was also temperature dependent. These 
observations are in accord with what is known about 
the normal state of underdoped Y123 from other experi- 
ments. 

The properties of the sharp bosonic resonance in the 
infrared scattering rate agree with the properties of the 
(tt, 7r) spin-flip neutron resonance mode. The frequency 
of the bosonic resonance 248 cm" 1 or 31 meV agrees 
with the frequency of the neutron mode measured by 
Stock et al. [2jj (provided the pseudogap is included in 
the model) . The temperature dependence of the strength 
of the mode also agrees with the neutron measurements 
on Ortho II Y123 samples from the same source. The 
total spectral weight of the mode determined from least 
squares fits drops monotonically with temperature and 
reaches zero between 200 and 300 K, in agreement with 
the spectral weight of the neutron mode. 

Fig El shows separately the two contributions to 
the scattering rate. We see that the combination of 
a sharp mode whose strength decreases with tempera- 
ture and a temperature independent background solves 
a long standing problem in the optical conductivity of 
the cuprates: there is a strong temperature dependence 
at low frequencies but above 500 cm -1 the scattering is 
nearly temperature independent. In the picture of mode 
and background, both caused by coupling to bosons, this 
high frequency behavior is the result of a competition be- 
tween a rising absorption due to occupancy of the back- 
ground boson population and a decrease in the strength 
of the bosonic mode with temperature. This scenario 
predicts an overall linear temperature dependence at all 
frequencies above the temperature where the sharp mode 
disappears. There is some evidence for this kind of be- 
havior in the data for La2- a; Sr a ;Cu04 which does not have 
a magnetic resonance but shows a linear temperature de- 
pendence characteristic of bosons both in the dc resistiv- 
ity [H3 and optical scattering rate at high frequency . 



FIG. 16: Temperature dependence of the amplitude of the 
sharp mode. The open diamonds with cross are from the 
second derivative analysis of the scattering rate, the closed 
hexagons are from fit a scattering rate to the model including 
a sharp mode and a background. The upright triangles show 
the energy integrated amplitude of the neutron mode Stock 
et al. |27| . 



Summary and conclusion 

In summary, we have shown that two optical phenom- 
ena in Ortho II Y123 have a close association with the 
magnetic resonance as measured for samples from the 
same source. The first is the onset of scattering associ- 
ated with a bosonic resonance. This has been suggested 
by several previous investigators 0,0,O,H3, but here, 
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for the first time we base our conclusions on experiments 
that use crystals with the same, well-defined doping level 
grown by the same technique. 

The second is the low-temperature conductivity peak 
that we have designated SI. Its connection to the mag- 
netic resonance is not obvious since the magnetic res- 
onance by itself is not optically active. One possible 
mechanism is the growth of in-plane coherence that takes 
place in this temperature region and has been suggested 
to be the closely related to the transverse plasma reso- 
nance that is seen in interplane conductivity |69j . The 
peak SI may be such a resonance associated with in- 
plane char ge d ensity modulations either through pinning 
to defects [lflj or through a transverse plasma resonance 
associated with the inhomogeneous charge density that 
seems to be p resent in underdoped cuprates in the normal 
state HSiJ. 
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